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Pseudoscalar Higgs boson production at hadron colliders in NNLO QCD 
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We compute the total cross-section for direct production of tiie pseudoscalar Higgs boson in 
hadron collisions at next-to-next-to-leading order (NNLO) in perturbative QCD. The C»(q^) QCD 
corrections increase the NLO production cross-section by approximately 20 — 30 per cent. 
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I. INTRODUCTION 

Supersymmetry is one of the most popular exten- 
sions of the Standard Model (SM). Generically, su- 
persymmetric theories predict a very rich spectrum of 
elementary particles. In particular, the Higgs boson 
sector of the Minimal Supersymetric Standard Model 
(MSSM) consists of two complex Higgs doublets. Af- 
ter electroweak symmetry breaking, three Higgs fields 
are absorbed by the and Z bosons into their longi- 
tudinal degrees of freedom; the remaining five degrees 
of freedom are physical Higgs bosons. In addition to 
the standard Higgs boson (/i) , a heavier neutral Higgs 
boson (H), two charged scalar Higgs bosons (H^), 
and a neutral pseudoscalar Higgs boson (A) appear in 
the spectrum. 

The tree-level masses of the Higgs bosons in the 
MSSM are usually described in terms of two inde- 
pendent parameters: the mass of the pseudoscalar 
Higgs boson tua and the ratio of the vacuum ex- 
pectation values of the two Higgs doublets tan/3 = 
vi/v2. Currently, these parameters are restricted by 
the experiments at LEP which set a lower bound 
niA > 91.9 GeV and exclude the values 0.5 < tan/3 < 
2.4 [1, 2]. Future searches for the pseudoscalar Higgs 
boson will be carried out at the Tevatron and at the 
LHC. It is therefore important to obtain a reliable 
theoretical estimate of its production cross-section in 
hadron collisions. 

The pseudoscalar Higgs boson does not couple to 
the gauge bosons at tree level; the major mechanisms 
for producing it are the gluon-gluon fusion mecha- 
nism gg ~f A and the associated production process 
99 J 99 ^99- The relative significance of the produc- 
tion channels depends on the mass of the Higgs boson 
and the value of tan/3. For larger values of tan/3, 
the coupling of the pseudoscalar Higgs boson to up 
quarks, ^up w„/ tan/3, decreases while the coupling 
of the Higgs boson to down quarks, gdown ~ rnd tan /3, 
increases. As a consequence, the phenomenology of 
the axial Higgs for large values of tan /3 is very differ- 
ent from the phenomenology of the SM Higgs boson. 



For the {t, b) family, the Higgs boson interaction with 
bottom quarks dominates for values of tan/3 > 10. 

Gluon-gluon fusion through a quark triangle-loop 
is the dominant production channel of the light pseu- 
doscalar Higgs boson in hadron collisions; there is no 
squark-loop contribution to the ggA coupling. In this 
Letter we study the gluon-gluon fusion cross-section 
for small and moderate values of tan /3. We can then 
neglect the contribution of bottom quarks and focus 
on the production of the pseudoscalar Higgs boson 
through the top-quark loop. In addition, we consider 
small values of the pseudoscalar Higgs boson mass, 
ruA < 300 GeV. Since niA ^ 2mt, where mj is the 
mass of the top-quark, the interaction of the pseu- 
doscalar Higgs boson with gluons and light quarks can 
be described by the effective Lagrangian [3] 
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The Wilson coefficients Gi and G2 are given by 
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where as(/^) is the MS coupling constant, Uf is 
the number of massless quark flavors and Lt = 
\og{ti^/mj). 

The renormalization in higher orders of the pertur- 
bation theory of the effective Lagrangian in Eq. (1) is 
subtle. Because of the axial anomaly, the derivative 
of the axial current of light quarks, d^J^ , mixes un- 
der renormalization with the operator G^i, G^'^ . The 
renormalization procedure must preserve the anomaly 
relation 
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A detailed discussion of the effective Lagrangian and 
its renormalization can be found in Ref. [3]. 

The Levi-Civita tensor e^i,ap and the 75 Dirac ma- 
trix are four- dimensional and their treatment in di- 
mensional regularization is a delicate problem. For 
calculational convenience we use the approach sug- 
gested in Ref. [4]. According to this prescription the 
75 matrix is represented as 75 = ■j/24 e^u^apl'^l'^j"''!^^ 
and the axial- vector current as = l/2'0(7^75 — 
With this substitution we can factor out the 
product of two Levi-Civita tensors from the produc- 
tion cross-section and evaluate it in terms of d = 4 — 2e 
dimensional metric tensors (5^jy, using 
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We have checked that the above prescription is con- 
sistent with the renormalization procedure of Ref. [3] 

by calculating the decay rate of the pseudoscalar 
Higgs boson through NNLO. Our results are in agree- 
ment with the expressions for the decay rate given in 
Ref. [3], where a four-dimensional treatment of the 
Levi-Civita tensors was employed. 

In this Letter we present the NNLO QCD cor- 
rections to the pseudoscalar Higgs boson production 
cross-section in hadron collisions. Various partonic 
processes contribute to the cross-section at this order. 
Specifically, we have to compute: a) virtual correc- 
tions to gg A, qq ^ A up to 0(al); b) virtual 
corrections to single real emission processes gg Ag, 
qg M, qg M, QQ Ag, up to 0(as); and c) 
double real emission processes gg Agg, gg Aqq, 
qg Aqg, qq Agg, qq Aqq. We evaluate 
the above corrections using the method introduced in 
Ref. [5] for the algorithmic evaluation of phase-space 
integrals. 



II. PARTONIC CROSS-SECTIONS 



Higgs boson [5] . We can then write 
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For convenience, we present here the difference 
^'Pij'' = ^ij^ ^ Vij^- The rjlj^ terms are listed in Sec- 
tion IV of Ref. [5]. 

We set the renormalization and the factorization 
scales equal to the mass of the Higgs boson niA- At 
LO we find 



S<l^ = 0. 



At NLO we obtain 
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The NLO terms are in agreement with the results of 
Ref. [6]. 

(2) 

At NNLO the differences S(plj are very simple. For 



both and ry-^- we find the same polylogarithmic 

terms; the differences contain simple logarithms and 
the Riemann ^2 constant. We obtain 
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In this section we present analytic expressions for 
the partonic cross-sections A+X, where i,j = 

q, q, g. We write 



an = a. 
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The coefScients ^-^^ are very similar to the coeffi- 
cients "rj^^^ of the perturbative expansion of the par- 
tonic cross-sections for the production of the scalar 
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The above results are valid if the renormalization and 
factorization scales are equal to the mass of the Higgs 
boson, /If = Hr = niA- The complete functional de- 
pendence of the partonic cross-sections on these scales 
can be easily restored by solving the DGLAP equation 
and the renormalization group equation and using the 
above expressions as the boundary conditions. This 
procedure is outlined in Ref. [5]. 



III. NUMERICAL RESULTS 
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FIG. 1: The pseudoscalax Higgs boson production cross- 
section at the LHC at leading (dotted), next-to- leading 
(dashcd-dottcd) and ncxt-to-ncxt-to-lcading (solid) order. 
The two curves for each case correspond to /Xr = M/ = 
m>i/2 (upper) and fj,r = fJ-f = 2mA (lower). 

We now discuss the numerical impact of the NNLO 
corrections on the pseudoscalar Higgs boson produc- 
tion cross-section at the LHC and the Tcvatron. To 
calculate the cross-section we must convolute the hard 
scattering partonic cross-sections of Section II with 
the appropriate parton distribution functions, accord- 
ing to the factorization formula 
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Here fl^^ is the MS distribution function of the parton 
i in the hadron h, ® denotes the standard convolution, 
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and X = m\/s, where s is the square of the total 
center of mass energy of the hadron-hadron collision. 

The complete NNLO parton distribution functions 
are not yet available. In Ref. [8] an approximate 
NNLO evolution [9] has been implemented in order 
to determine the NNLO MRST parton distribution 
functions. We use these approximate solutions for the 
numerical evaluation of the pseudoscalar Higgs boson 
production cross-section keeping the same initializa- 
tion parameters as in Ref. [5]. 

To demonstrate the convergence properties of the 
perturbative series for the hadronic cross-section, we 
present the LO, NLO and NNLO results for both the 
LHC and the Tevatron. In order to improve upon the 
heavy-top quark approximation, we normalize our re- 
sults to the leading-order cross-section with the exact 
dependence on m^. 

The total cross-section for the LHC is shown in 
Fig. 1. From Fig.l we observe that the scale de- 
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FIG. 2: The pseudoscalar Higgs boson production cross- 
section at the Tevatron at leading (soUd), next-to- leading 
(dashed) and next-to-next-to-leading (dotted) order. The 
two curves for each case correspond to jir = = vhaI^ 
(upper) and fir = Hf = 2mA (lower). 

pcndence of the Higgs production cross-section at 
NNLO is approximately 15%; this is a factor of two 
smaller than the NLO scale dependence and a fac- 
tor of four less than the scale variation at LO. Despite 
the scale stabilization, the corrections are rather large; 
the NLO corrections increase the LO cross-section by 
about 70%, and the NNLO corrections further in- 
crease it by approximately 30%. The K factor, de- 
fined as the ratio of the NNLO cross-section and the 
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LO cross-section, is approximately two. In Fig. 2 we 
plot the values of the Higgs production cross-section 
at the Tevatron. The NNLO K factor is approxi- 
mately three, and the residual scale dependence is 
approximately 30%. The K factors for the produc- 
tion of the pseudoscalar Higgs boson and the scalar 
Higgs boson [5, 7] are comparable in magnitude; this 
is a consequence of the similarity of the corresponding 
partonic cross-sections as discussed in Section 11. 
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FIG. 3: The pseudoscalar Higgs boson production cross- 
section at the LEG at loading (dotted), next-to-leading 
(dashed) and next-to-next-to-leading (solid) order as the 
function of factorization and renormalization scale /tt. The 
mass of the Higgs boson is 115 GeV. 

In Ref. [5] we have argued that since the dominant 
contribution to the integrated cross-section for the 
scalar Higgs boson comes from the region close to 
the Higgs boson production threshold, we should 
choose values of the scale n which are smaller 
than the mass of the Higgs boson. This choice 
decreases the NNLO corrections and the Higgs boson 
production cross-section increases as compared to 
conventional choice of the scales, = l-^f = nT-H- 
The production cross-section for the pseudoscalar 
Higgs boson exhibits the same behavior. In Fig. 3 
we show an example, whore we plot the production 
cross-section for niA = 115 GeV. We equate the 
renormalization and factorization scales and vary the 
factorization scale from ^ = 25 GeV up to the mass 
of the pseudoscalar Higgs boson. The plot illustrates 
that for smaller values of /U, the NLO cross-section 



increases more rapidly than the NNLO cross-section, 
and the difference between the NLO and the NNLO 
results becomes smaller. Therefore, the convergence 
of the perturbative series is improved for smaller 
values of the factorization scale. 

IV. CONCLUSIONS 

In this Letter we presented the NNLO QCD cor- 
rections for the production cross-section of the pseu- 
doscalar Higgs boson in hadron collisions. Our results 
are valid in the heavy top-quark limit and for small 
to moderate values of tan/3. 

The analytic expressions which we presented here 
and those for the scalar Higgs boson [5] are very simi- 
lar. In both cases, the QCD corrections are large and 
important for quantitative estimates of the hadronic 
cross-sections at the Tevatron and the LHC. The size 
of the NNLO corrections for the pseudoscalar Higgs 
boson indicates that the perturbative expansion of the 
production cross-section converges, albeit slowly. A 
similar convergence behavior was observed for the SM 
Higgs boson hadronic production cross-section [5, 7]. 

In order to verify the compatibility of the prescrip- 
tion of Ref. [4] for the Levi-Civita tensor with the Wil- 
son coefficients for the effective Lagrangian in Eq. (1), 
derived in Ref. [3], we computed the pseudoscalar 
Higgs boson decay rate through NNLO in QCD. Our 
results are in complete agreement with the expressions 
for the decay rate in [3] . 

We note that the calculation reported in this Let- 
ter has been performed using the nic^thod of Ref. [5] 
which combines the optical theorem with integration- 
by-parts reduction algorithms to achieve a systematic 
evaluation of phase-space integrals. As our calcula- 
tion demonstrates, the method is general and process- 
independent. We are confident that the same method 
will be very useful in studying other processes of phe- 
nomenological interest. 

Finally, as we completed this manuscript, we be- 
came aware of a similar calculation by R. Harlander 
and W. Kilgore [10]. We have compared results for the 
partonic cross-sections and find complete agreement. 
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